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Abstract 

We extend the semiclassical theory of scarring of quantum eigenfunctions ip n (q) 
by classical periodic orbits to include situations where these orbits undergo generic 
bifurcations. It is shown that IVVilq)) 2 , averaged locally with respect to position q 
and the energy spectrum {E n }, has structure around bifurcating periodic orbits with 
an amplitude and length-scale whose /i-dependence is determined by the bifurcation 
in question. Specifically, the amplitude scales as h a and the length-scale as W , 
and values of the scar exponents, a and u, are computed for a variety of generic 
bifurcations. In each case, the scars are semiclassically wider than those associated 
with isolated and unstable periodic orbits; moreover, their amplitude is at least as 
large, and in most cases larger. In this sense, bifurcations may be said to give rise to 
superscars. The competition between the contributions from different bifurcations 
to determine the moments of the averaged eigenfunction amplitude is analysed. We 
argue that there is a resulting universal ^-scaling in the semiclassical asymptotics of 
these moments for irregular states in systems with a mixed phase-space dynamics. 
Finally, a number of these predictions are illustrated by numerical computations for 
a family of perturbed cat maps. 



* Short title: Bifurcations and scarring 
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1 Introduction 



One of the main goals in quantum chaology has been to determine the link between 
classical periodic orbits and quantum spectral fluctuations in the semiclassical limit. In 
fully chaotic systems, where the periodic orbits are isolated and unstable, this connection 
is embodied in Gutzwiller's trace formula (Gutzwiller 1971), and in integrable systems by 
a corresponding expression involving the phase-space tori (Berry & Tabor 1976). These 
formulae fail, by diverging, when periodic orbits bifurcate; that is, when combinations of 
stable and unstable orbits collide and transmute, or annihilate, as a system parameter 
varies - phenomena that characterize dynamics in the mixed regime. They must then 
be replaced by transitional or uniform approximations which interpolate through the 
bifurcation (Ozorio de Almeida & Hannay 1987; Tomsovic et al. 1995; Ullmo et al. 1996; 
Sieber 1996; Schomerus & Sieber 1997; Schomerus 1998; Sieber & Schomerus 1998). 

That individual orbit bifurcations can have an important, and sometimes dominant 
influence on spectral statistics was pointed out by Berry et al. (1998), and demonstrated 
for a particular example, the perturbed cat maps. More generally, Berry et al. (2000) 
developed a semiclassical theory for the competition between the various generic bifur- 
cations found in Hamiltonian systems to determine the moments of the quantum energy 
level counting function. This suggests that these moments diverge in a universal way, 
characterized by certain twinkling exponents, as h — > 0. 

A second major goal of quantum chaology has been to understand the influence of 
classical periodic orbits on quantum wavefunctions in the semiclassical limit. It was first 
noticed by McDonald (McDonald 1983) that individual eigenfunctions can have enhanced 
intensity along short periodic orbits in classically chaotic systems. This phenomenon was 
later studied systematically by Heller (Heller 1984), who called such structures scars. He 
developed a theory of scarring, based on wavepacket dynamics, which has subsequently 
been extended to describe a variety of statistical properties of quantum chaotic eigenfunc- 
tions (Kaplan 1999). 

An alternative theory of scarring, based on an approach closely related to the trace 
formula, was initiated by Bogomolny (Bogomolny 1988). In this, the semiclassical ap- 
proximation to the energy-dependent Green function is used to show that for quantum 
eigenfunctions ^ n (q) corresponding to energy levels E n in a fixed energy range, (|^ n (q)| 2 ), 
where (• • • ) denotes an average over the states in question and locally over position q, has 
complex-Gaussian fringes with, in two-degree-of-freedom systems, amplitude and length- 
scale of the order of h 1 ^ 2 around unstable periodic orbits. A corresponding theory for 
Wigner functions was developed by Berry (1989). 

We emphasize two limitations of the theories mentioned above. First, they only de- 
scribe scarring in eigenfunctions that have been averaged over an energy interval which, 
semiclassically, contains a large number of states. Resummation techniques have been 
applied to provide some information about individual eigenfunctions (Agam & Fishman 
1994; Fishman et al. 1996), but a detailed understanding of the phenomenon in this 
case remains to be developed. Second, they concentrate specifically on the influence of 
periodic orbits which are isolated and unstable. (Semiclassical theories describing the 
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connection between quantum wavefunctions and phase-space tori in classically integrable 
systems have also been developed; see, for example, Berry 1983 for a detailed review.) 

Our purpose here is to address the second of these limitations. Specifically, our aim 
is to show how Bogomolny's theory can be extended to include the description of semi- 
classical structures in quantum eigenfunctions associated with generic classical periodic 
orbit bifurcations in systems with two degrees of freedom. We focus in particular on 
the ^-dependence of the amplitude and length-scale of the fringes corresponding to those 
identified by Bogomolny. Our main result is that the amplitude is of the order of h a , 
and the length-scale is of the order of H u , where a and u are bifurcation-dependent scar 
exponents whose values we calculate in a number of different cases. Crucially, u < 1/2 
for all the bifurcations studied, and for most a < 1/2 as well. In this sense, bifurcations 
may be said to give rise to superscars. In order to quantify this, we determine the way in 
which bifurcating orbits contribute, via a competition, to the semiclassical asymptotics 
of the moments of (|^n(<?)| 2 ), in the same way as was done for spectral fluctuations by 
Berry et al. (2000). It is argued that this competition results in universal ft-scalings of the 
moments for the irregular eigenfunctions of systems with mixed phase-space dynamics. 
Finally, as an example, we apply some of the techniques developed to study the influence 
of one particular bifurcation on the eigenfunctions of a family of quantum perturbed cat 
maps. 



2 Scar Formulae 

Our aim in this section is to derive semiclassical scar formulae for bifurcating periodic 
orbits which generalize those obtained in Bogomolny (1988) for unstable orbits far from 
bifurcation. 

We begin, following Bogomolny, with the energy-dependent Green function 

G(q , q;E) = E «MM 3 ) i (1) 

n 

where ip n (c\) is the eigenfunction of the quantum Hamiltonian corresponding to the energy 
level E n . The identity we seek to exploit follows from setting q' = q: 

V |^ n (q)| 2 5 £ (E - E n ) = -- ImG(q, q; E + ie). (2) 
z — <■ n 

n 

Here, 5 e (x) is a normalized, Lorentzian-smoothed 5-function of width e. (It is straightfor- 
ward to transform (pj) to give differently smoothed ^-functions, for example Gaussians.) 
The left-hand side of (pj) thus corresponds to a sum over eigenstates for which E n lies 
within a range of size of the order of e centred on E. Semiclassically, it is approximately 
the average of |^ n (q)| 2 with respect to these states multiplied by d(E), the mean level 
density. For systems with two-degrees-of-freedom 

3(b » ~ S <3) 
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as h — > 0, where 

V(E) = J 5{E- H(p, q))d 2 qd 2 p (4) 

and H (p, q) is the classical Hamiltonian. 

The connection with classical mechanics is achieved using the semiclassical approxi- 
mation to the Green function. For systems with two-degrees-of-freedom, this is 

G(q', q; E) « T^= E v^ eXp {^ q '' q; E ^ ~ f ( 5 ) 

where the sum includes all classical trajectories from q to q' at energy E, S" 7 is the action 
along the trajectory labelled 7, 

a 2 s^ a 2 s~, 



D, = det[ |gf , (6) 

\ asaq as 2 / 

and v is the Maslov index (Gutzwiller 1990). When q' = q, the sum in (0) is clearly over 
closed orbits. 

We note in passing that it follows from J2J that 

^5 e (E - E n ) = —1m ! G(q,q;E + ie)d 2 q. (7) 

n J 

Substituting in the closed orbit sum for G(q, q; E+ie) and integrating term-by-term using 
the method of stationary phase leaves contributions from the periodic orbits. Assuming 
these are all isolated, as is the case for hyperbolic systems, the result is the trace formula 
(Gutzwiller 1971) 



p r=l 



T 

X 



det (M;-i)| 



cos ( - ru p - ] <>xp 



rT p 



(8) 



where p labels primitive periodic orbits with period T p and monodromy matrix M p , and r 
labels repetitions. As noted in the Introduction, this formula fails at bifurcations, where 
det (Mp — I) = 0. Assuming that the periodic orbits lie in families which form tori in 
phase space gives the corresponding expression for integrable systems (Berry & Tabor 
1976). 

Bogomolny's scar formula follows not from integrating over all positions q, as in (J2J), 
but from performing a local average of ((2|) with respect to q (we postpone specifying the 
size of the averaging range until after the result has been stated), which we take to be 
smooth (convolution with a normalized Gaussian, for example). On the left-hand side this 
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gives, approximately, d(E) (|^ n (q)| 2 ), where (• ■ ■ ) denotes a combination of the spectral 
average described above and the local q-average. On the right-hand side, the q-average 
selects from the closed orbits those that are close to periodic orbits (i.e. for which the 
change in momentum after return is appropriately small). These can then be described 
by linearizing about the periodic orbits. Essentially, this corresponds to expanding the 
action up to terms which are quadratic in the distance from the periodic orbit. The result 
is that 



n 

1 



n(q;£)- 



7T 




r=l \Z\ 



X 



[M;( 



[2 



cxp 



l det(M;-l) 

M rtp+ 2 [m;(,)] 12 y 



7T 



cxp 



— e- 



rTp 

h 



(9) 



where z is a coordinate along a given periodic orbit and y is a coordinate transverse to 
it, [MJ(z)]^. denotes the elements of the monodromy matrix (which are functions of z), 
z is the velocity along the periodic orbit, and 



fi(q;£)= / 5(E-H(p,q))d 2 p. 



(10) 



This in turn implies that 



>n(q)| 2 > 



Q(q;E) iTrVh 
V(E) ~V(E) 



Im 



exp 



, / ldet(M;-l) 
h V bp+ 2 




cxp 



— e- 



rTp 
h 



(11) 



Equation (fTT| is Bogomolny's scar formula. In classically ergodic systems, the first 
term represents the quantum-ergodic limit of the eigenfunction probability density (Shnirelman 
1974, Colin de Verdiere 1985, Zelditch 1987). The second describes complex Gaussian 
fringes (the y-dependent part), with length-scale and amplitude both of the order of ft 1//2 , 
associated with each periodic orbit. This structure will be resolved if the local q-average 
is over regions whose dimensions are small compared to the length-scale of the fringes; 
that is, over regions whose dimensions scale as h 6 , where 5 > 1/2. In order for the near-to- 
periodic-orbit approximation to be valid, we must also have 5 < 1; that is, the dimensions 
of the averaging range must be large compared to a de Broglie wavelength. 

The trace formula © can be recovered from Q by integrating over z and y. The 
z-integral gives the period in the amplitude of the periodic orbit contributions, and the 
?/-integral gives the determinant. Note that the power of H in the trace formula amplitude 
is the amplitude exponent in (JUJ), —3/2 (which in turn is equal to the amplitude exponent 
in (|TT| minus two - the exponent in d), plus the length-scale exponent of the fringes, 1/2. 
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The approximations © and (fTTj) break down in two ways. First, at self-focal points 
along an orbit [M£( z)\ 12 = and the amplitude diverges. This can be remedied straight- 
forwardly using Maslov's method, and we will not concern ourselves further with it here. 
Second, when an orbit bifurcates det (M£ — I) = 0, and so the formulae become y- 
independent. Essentially, this means that the fringes are infinitely wide (it is this infinity 
which, upon integration with respect to y, transfers itself to the amplitude in the trace 
formula). Our purpose in this paper is to show how to correct (j^J) and (fTTj) in this case. 

It might be thought that the scar formulae for bifurcating orbits could be obtained eas- 
ily by expanding the action in (jSJ) to higher order than quadratic. For some of the simpler 
bifurcations (e.g. the codimension-one bifurcations of orbits with r = 1 and r = 2) this is 
correct (see the example in Section 4). However, for more complicated bifurcations it is 
incorrect, because for these the linearized map is equal to the identity, which cannot 
be generated by the action S(q', q; E). Thus it is difficult to build into the semiclassical 
expression for the qq'-representation of the Green function a well-behaved description 
of the nonlinear dynamics which the linearized map approximates. The solution to this 
problem, originally proposed in Ozorio de Almeida & Hannay (1987), is to transform the 
Green function to a mixed position-momentum representation, and this is the approach 
we now take. 

The Ozorio de Almeida-Hannay method involves, first, Fourier transforming G(q', q; E) 
with respect to q'. This gives the Green function in the q p'-representation, G(p', q; E) 
(p ; is the momentum conjugate to q'). The semiclassical approximation to G takes the 
same form as (pj), except that S(q', q; E) is replaced by the q p'-generating function 
S(p',q;E). G(q', q; E) may then be rewritten, semiclassically, as the Fourier transform 
of this expression with respect to p'. (For an alternative approach leading to the same 
final answer, see Sieber 1996). The result, for the semiclassical contribution to G(q, q; E) 
from closed orbits in the neighbourhood of a bifurcating periodic orbit, takes the following 
form. 

Consider the case of a codimension-K bifurcation of a periodic orbit with repetition 
number r. As before, let z be a coordinate along the orbit at bifurcation, let y be 
a coordinate transverse to it, and let p y be the momentum conjugate to y, so that y 
and p y are local surface of section coordinates. Let <& rt ic(y,py,it) be the normal form 
which corresponds to the local (reduced) generating function in the neighbourhood of 
the bifurcation (Arnold 1978; Ozorio de Almeida 1988), where x = (xi,x 2 ,--- , Xjc) are 
parameters controlling the unfolding of the bifurcation. Then, up to irrelevant factors, 
the contribution to G(q, q; E) is 



G r>K (y;x) = J exp 



dp y (12) 



(as already stated, we are here interested in determining the ^-dependence of the ampli- 
tude and length-scale of the associated fringes, and so have neglected terms in (fT2j) . such 
as an ^-independent factor in the integrand, which do not influence these). 

Before proceeding further, we make three remarks about (fT2|) . First, the power of h 
outside the exponential arises from adding 1/2, which comes from the Fourier transform, 
to the exponent in 3/2. Second, the representation used in Berry et al. (2000) for 
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the fluctuating part of the spectral counting function may be derived from (fT2|) by taking 
the trace of G r> K, which involves integrating the right-hand side of (fT2|) with respect to y 
(the z-integral is trivial, as before), and then integrating with respect to E resulting in 
a further multiplication by h. Likewise, the formulae of Ozorio de Almeida and Hannay 
(1987) correspond to taking the trace of G t% k, keeping the terms we have neglected. Third, 
the ^-dependence of the fringes in Bogomolny's scar formula for unstable periodic orbits 
far from bifurcation can be recovered using the appropriate normal form: 

$r,0 =P 2 y + V 2 , (13) 

which corresponds to a particular, ^-independent choice of units for p y and y. Evaluating 
the integral then gives 

G>,o(q) oc exp , (14) 



as in p. 

Equation (fT2|) is the starting point for the analysis of bifurcating orbits. Our strategy 
is essentially the same as that used in Berry et al. (2000) to study the related fluctuations 
in the spectral counting function (see also Berry 2000 for a review of applications to other 
areas in wave physics): first, rescale y and p y to remove the 1/h factor from the dominant 
term (germ) of Q r ,K in the exponent, and then apply a compensating rescaling of the 
parameters x±, X2, ■ • ■ , xk to remove the ^-dependence from the other terms which do not 
vanish as h — > 0. This will lead to 

G r , K (y; x; h) = j^G^ K (jJL., {^^} , 1) • (15) 

The exponent a describes the semiclassical amplitude of the fringes in (|?/> n (q)| 2 ) as- 
sociated with the bifurcation, and the exponent ui describes the ^-dependence of their 
length-scale, or width. We call these the scar exponents. Note that the corresponding 
amplitude exponent in J2 n (l^ni^l 2 ) q 5 E {E — E n ) is 2 — a. The exponents a describe the 
range of influence of the bifurcation in the different unfolding directions x n . Their sum 

K 



lr,K — Vn,r,K (16) 



n=l 



describes the ^-scaling of the ^-dimensional x-space hypervolume affected by the bifur- 
cation. 

We now calculate these exponents in a variety of examples. Consider first the r = 1 
bifurcations which correspond to cuspoid (i.e. corank 1) catastrophes. For these, the 
normal forms are (Berry et al. 2000) 

K 

$ ltK (y,p y ;x) =p 2 y + y K+2 + ^x n y n . (17) 

n=l 



7 



Substituting this into (fT2|) and evaluating the integral then gives 



G 1>K (y;x) oc ^ ex P 




K 



K+2 



+ XnyT 



n=l 



(18) 



Making the rescalings y = y/h l ^ K+2 \ x n = x n /h l n /( K+2 ) 
the exponent, and so for the cuspoids we have 



removes the ^-dependence of 



1 

2' 



K + 2' 



n 



K + l 



(19) 



and 



ll,K 



K(K + 3) 
2{K + 2) - 



(20) 



Analogous expressions can be written down for the r = 1 bifurcations corresponding to 
the more complicated case of catastrophes of corank 2 (see, for example, Berry 2000) in 
the same way. 

When r > 1, the generic bifurcations with K = 1 have been classified by Meyer (Meyer 
1970, 1986; Arnold 1978; Ozorio de Almeida 1988), and those with K = 2 by Schomerus 
(1998). The relevant parts of the corresponding normal forms, taken from Berry (2000) 
(to which readers are referred for further details), are summarized in Table 1 (in the 
expressions given, we are retaining only those terms which affect the exponents we seek 
to calculate). 



r,2 



2 
3 
4 

5 
> 6 



(p 2 y + y 2 ? 
p 2 y y 2 



+ x 1 (p, 
4- xi (p 2 



v + V 



Re [{p y + iy) 

ipl + y 



2\3 



y . 

+ x\{p 2 

2 



) + x 2 Re [{p y + iy)' 

+ x 2 (p 2 v -y 2 ) 2 



y + y 2 )+x 2 (p 2 y 
+ x 1 (p 2 y + y 2 ) +x 2 (p 2 y + y 



y 

2\2 



2\2 



Table 1: The relevant parts of the normal forms for K = 2 bifurcations of period-r orbits 
(taken from Berry 2000). The corresponding expressions for K = 1 bifurcations, & r ,i, 
follow from settling x 2 = 1. 

The corresponding scar exponents, and the hypervolume exponents 7 are given in 
Table 2 (K = 1) and Table 3 (K = 2). 

Finally, we consider bifurcations of orbits for which r > 2K + 2. In this case, the 
relevant terms in the normal forms are (Berry et al. 2000) 

K 

$,,*(</, p,;x) = I K+l + J2*nI n + 0(I K+2 ), (21) 

n=l 
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r 


Oi r i 




7r 1 


2 


1/2 


1/4 


1/2 


3 


1/3 


1/3 


1/3 


> 4 


1/4 


1/4 


1/2 











Table 2: Scar exponents for generic, codimension-1 bifurcations. 



r 


a r ,2 


UV,2 


7r,2 


2 


1/2 


1/6 


1 


3 


1/4 


1/4 


3/4 


4 


1/4 


1/4 


1/2 


5 


1/5 


1/5 


4/5 


> 6 


1/6 


1/6 


1 



Table 3: Scar exponents for generic, codimension-2 bifurcations. 

where 

i = y 2 + vl- (22) 

Expressing <3> in terms of y and p y , we find 

1 

®r,K - ^r,K - - ^ (23) 

and 

lr,K = y (24) 

(c.f. the r > 4 exponents in Table 2, and the r > 6 exponents in Table 3). 

The main point we wish to draw attention to is that, in all the cases listed above, 
uj < 1/2 and a < 1/2, and that in most cases a < 1/2. Recall that ui = a = 1/2 for 
periodic orbits far from bifurcation. In this sense, bifurcations may be said to give rise to 
superscars; that is, to scars that are semiclassically wider, and often greater in amplitude 
than those associated with non-bifurcating orbits. We shall demonstrate this with an 
explicit example in Section 4. 

The width exponent u determines the scale for the q-average in (|^ n (q)| 2 ) which 
allows the fringe structure to be resolved. Specifically, if the dimensions of the averaging 
range scale semiclassically as h 5 , the fringe structure will be resolved if 5 > u, but not 
if 5 < u. Recall that for the approximation to hold in which the main contribution to 
the average comes from closed orbits in the neighbourood of periodic orbits, we must also 
have 5 < 1; that is, the average must extend over many de Broglie wavelengths. 

We also note that the scar exponents satisfy 

P = l-a-u, (25) 
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where (3 is the amplitude exponent of the fluctuations in the spectral counting function 
associated with the bifurcation in question (Berry et al. 2000). This follows from a 
comparison of (fT2|) with the corresponding expression for the counting function, which, as 
already noted, corresponds to integrating (fT2l) with respect to y and multiplying by h. It 
generalizes the connection discussed above between the power of h in the trace formula for 
d(E) and the amplitude and width exponents of Bogomolny's fringes for non-bifurcating 
orbits. 

3 Moment asymptotics 

One way to quantify scarring effects is in terms of the moments of the wavefunctions. 
Consider the case when all periodic orbits are isolated and unstable. We define 

c 2m (h) = ±J ((|^(q')l 2 ) - ^p) 2 "^V (26) 

where the q'-integral is over an ^-independent volume Aq. Note that these are the mo- 
ments not of the amplitude of the wavefunction itself, but of the amplitude averaged 
with respect to position (over a region which shrinks as h — > 0, but which contains an 
increasing number of de Broglie wavelengths) and energy (semiclassically many levels) . 

Assuming that the wavefunctions are quantum ergodic on the scale of the local q- 
average implies that C 2m (ft) — > as K —* when m > 0. It follows from the fact that the 
fringes in (fTTf have scar exponents a = uj = 1/2 that their individual contributions to the 
moments scale as fo m+1 / 2 in this limit. The corresponding contribution from a bifurcating 
orbit is of the order of fr 2ma+UJ ? and so is semiclassically larger. For the bifurcations of 
periodic orbits with r = 1 corresponding to the cuspoid catastrophes we have that 

fim,i,K = Zma 1)K + Ui tK = m+ 1 (27) 

A + A 

The values of these exponents for the generic bifurcations with K = 1 and K = 2 when 
r > 1 are given in Table 4 and Table 5, for m < 3. 



r 


fil,r,l 


^1^,1 


P<2,r,l 








2 
3 
> 4 


5/4 
1 

3/4 


7/4 
4/3 
5/4 


9/4 
5/3 
5/4 


11/4 

6/3 

7/4 


13/4 

7/3 
7/4 


15/4 

8/3 
9/4 



Table 4: Values of ^ m ,r,i = 2ma rj i + u r ^ and /i mrl = 2m« r]1 + o; nl + 7 rjl for the 
codimension-1 scar exponents listed in Table 2. 

The moments defined by (j26| are implicitly functions of the system parameters x. 
Averaging them with respect to x produces an opportunity for a competition between 
the various generic bifurcations. The contribution of each bifurcation must be weighted 
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r 




Ml,r,2 




^2,r,2 


fis,r,2 




2 


7/6 


13/6 


13/6 


19/6 


19/6 


25/6 


3 


3/4 


3/2 


5/4 


2 


7/4 


5/2 


4 


3/4 


5/4 


5/4 


7/4 


7/4 


9/4 


5 


3/5 


7/5 


1 


9/5 


7/5 


11/5 


> 6 


1/2 


3/2 


5/6 


11/6 


7/6 


13/6 



Table 5: Values of /J, m ,r,2 — 2ma r]2 + io r ^ and \i mr2 
codimension-2 scar exponents listed in Table 3. 



2ma r2 + w r 2 + 7r2 for the 



by the associated hypervolume in x-space, and so scales as fi 2ma + u +~i ; provided that the 
x-average of the 2mth power of the ^-independent term in (fl~5|) exists (see Berry 1977 for 
a discussion of this subtle point). For the r = 1 bifurcations corresponding to the cuspoid 
catastrophes, 

K + l , . 

/V,1,K = 2ma l,K + Ul,K + ll,K = TTl H . (28) 

When r > 1, the values of these exponents for the generic bifurcations with K = 1 and 
K = 2 are also listed in Table 4 and Table 5. The bifurcation that wins the competition, 
and hence which determines the rate at which the x-averaged moments tend to zero in 
the semiclassical limit, is the one for which 2ma + to + 7 is minimized; that is, 

1 (C 2m (h)) x = o(fr) (29) 



for any e < 0, with 



\i m = min(2ma + to + 7) 



(30) 



where (C2 m (h)) x denotes the x-averaged moments and the minimum is with respect to 
the generic bifurcations. This, of course, assumes that the minimum exists. We now argue 
that it does. 

Our reasoning is based directly on that of Berry et al. (2000), where the analogous 
problem of the moments of fluctuations in the level counting function was considered. 
First, we note that if the competition is restricted to bifurcations with r > 2K + 2, then 



min(2ma + uj + 7) = min^ 



2m + 1 



+ 



K 



(31) 



2{K + l) 

which exists for any m and can be calculated straightforwardly. Second, it was shown by 
Berry at al. that each bifurcation with r > 2K + 2 has a counterpart with r < 2K + 2 
with the property that the counterpart has a normal form with the same germ, and so the 
same a and uj exponents, but a larger 7 exponent. Hence (|3T|) represents the minimum 
with respect to all of the generic bifurcations and so 



2m + 1 K 
2(K+1) + ~2~ 



(32) 
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For example, fii = 5/4 (coming from K — 1), // 2 = 7/4 (also coming from K — 1), and 
/x 3 = 13/6 (if = 2). In general /i m sa y/2m + 1 - 1/2. 

It is natural to compare /i m to the corresponding exponent for Gaussian random 
functions, which are often taken as models of quantum chaotic wavefunctions. In that 
case, the moments (jZBj) are semiclassically of the order of h m . This follows from the 
results of Section IIIB of Eckhardt et al. 1995, if the operator considered there is the 
characteristic function of the region over which the local q-average in (|26|) extends. The 
same rate of vanishing also holds for the eigenvectors of random hermitian matrices (see 
Section IIIA of Eckhardt et al). (Readers are referred to Backer et al. 1998 for a detailed 
review of the rate of quantum ergodicity, and its characterization by moments analogous 
to those defined by (|26|) ). Crucially, we note that \i m < m for m > 1, and so then, if 
the background to the scars due to individual periodic orbits is modelled by a Gaussian 
random function, bifurcations dominate the semiclassical asymptotics. The contributions 
from individual non-bifurcating orbits are always subdominant. 

To summarize, the exponents /i m , which are analogous to the twinkling exponents of 
Berry et al (2000), determine the asymptotic scaling of the parameter-averaged moments 
C 2m in the limit as h — > when m > 1. Note that they are universal, that is, system 
independent. Note also that they are determined solely by generic bifurcation processes. 
As pointed out in the Introduction, these processes are characteristic of mixed phases- 
space dynamics, and so one might expect the exponents l|32*J) to describe the semiclassical 
deviations of the irregular (in the sense of Percival 1973) eigenfunctions in mixed systems 
from their ergodic limit. (They do not describe the regular eigenfunctions, for which the 
corresponding moments have a different origin, and can be calculated using the results of 
Berry et al. 1983.) 

4 Perturbed cat maps 

We now illustrate some of the general ideas described in the previous sections by focusing 
on a particular example: a family of perturbed cat maps. 
The maps we consider are of the form 

te)=0 ( :i ) + ^° s(2 ^> ) m ° d <33) 

where q and p are coordinates on the unit two-torus, and are taken to be a position and its 
conjugate momentum. These maps are Anosov systems for k < K max = (v3 — l)/v5 ~ 
0.33; that is, for k in this range they are completely hyperbolic and their orbits are 
conjugate to those of the map with k = (i.e. there are no bifurcations). Outside this 
range, bifurcations occur, stable islands are created, and the dynamics becomes mixed 
(see, for example, Berry et al. 1998, where these systems were used to demonstrate the 
influence of periodic orbit bifurcations on long-range spectral statistics). 

The quantization of maps like (|33f was developed by Hannay & Berry (1980), when 
k = 0, and Basilio de Matos & Ozorio de Almeida (1995) for non-zero k. The quantum 
kinematics associated with a phase space that has the topology of a two-torus restricts 
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Planck's constant to taking inverse integer values. The integer in question, N, is the di- 
mension of the Hilbert space of admissible wavefunctions. With doubly periodic boundary 
conditions (see, for example, Keating et al. 1999), these wavefunctions in their position 
representation have support at points q = Q/N, where Q takes integer values between 
1 and N. They may thus be represented by iV-vectors with complex components. The 
quantum dynamics is then generated by an iV x iV unitary matrix U whose action on the 
wavefunctions reduces to ijHHJ) in the classical limit; for example 



U, 



Qi,Q 2 



cxp 



2tti iN 

-rr(Q? - Q1Q2 + Ql) + —KsinpTTQx/N) 

iV Z7T 



(34) 



This matrix plays the role of the Green function of the time-dependent Shrodinger equa- 
tion for flows. 

Denoting the eigenvalues of U by e l9n , and the corresponding eigenfunctions by ^f n (Q), 
we have that 



J2\^n(Q)\ 2 S e (6-6 n ) = l + BB^2u^ Q exp(-i6k-ek) 



where 



k=l 



1 — e 6 cos x 



1 + e 2e — 2e e cos x 



(35) 



(36) 



is a periodized, Lorentzian-smoothed 5-function of width s (Keating 1991). Equation (I3H1) 
is the analogue for quantum maps of (j2J). The left-hand side corresponds, approximately, 
to N times the local n-average (over a range of size of the order of e) of |\l/ n ((5)| 2 , and so, 
dividing both sides by N and averaging smoothly with respect to Q (for example, taking 
the convolution with a normalized Gaussian) over a range large compared to a de Broglie 
wavelength (AQ = 1) but small compared to N, 



00 

*n(Q)| 2 ) « ^ + ^ Re ( U Q,q)q ex P (~ idk ~ £k )- 



(37) 



k=l 



Here (...) denotes a combination of the n-average and the Q-average (. . . )q. 

In our computations we took e large enough so that the dominant contributions to 
(j35|) and (|3*7j) come from the k — 1 terms in the sums on the right, and so, for example, 
we may substitute (l3~4"J) directly into (f37l) . In the semiclassical limit, as iV — > 00, the 
Q-average selects regions close to stationary points of the phase of (|3~31) . which we denote 
by Q/N = qf. These stationary points coincide with the positions of the fixed points of 
the classical map (|33|) : that is, qf satisfies 

1f = l(j-^^s(2nq f ^ (38) 

for integers j such that < g/ < 1 (see, for example, Boasman & Keating 1995). Ex- 
panding the phase of (j3*4"l) around qf gives, up to cubic terms, 



exp 



2mS f + ixiN(2 - Ksm(2nq f ))y 2 



2ir 2 iN 



k cos(2nq f)y 3 



(39) 
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where 



V = jj-Qf (4°) 

and 2nSf denotes the phase evaluated at g/. 

Provided that 2 — Ksin(27rgj) 7^ 0, this approximation is dominated by the quadratic 
term in the exponent when y is small. It thus describes complex-Gaussian fringes around 
the classical fixed points with a length-scale (in terms of y) of the order of N~ 1 / 2 . These 
are the analogues of Bogomolny's fringes. They will be resolved if the local Q-average is 
over a range that is small compared to N 1 / 2 (but which still grows as N — > 00). 

For the example being considered here, when k < k* ~ 5.94338 the two values of j 
in (JHHJ), j = and j = 1, each give rise to a single unstable fixed point for which the 
condition 2 — Ksin(27rg/) 7^ is satisfied. In Figure 1, we plot the left hand side of (|3~5|) 
when k = 3, with e = 2.2 and for N = 1597. The structure around the fixed points 
is clearly visible, and is most easily seen by applying the local Q-average (in this case, 
making a convolution with a normalized Gaussian of width 0.02). 

It is at bifurcations that 2 — Ksin(27rg/) = 0. Then the quadratic term in (|39|) vanishes, 
and the fringe structure comes instead from the cubic term. It thus has a ?/-length-scale 
of the order of iV~ 1//3 . The amplitude is the same as in the case of isolated fixed points 
{N~ 1 / 2 in the contribution to (\^ n (Q)\ 2 )). In the language of Section 2, this corresponds 
to a codimension-one bifurcation of a periodic orbit with r = 1 (a tangent bifurcation). 

In our example, the first bifurcation occurs when k = k*. At this parameter value, 
two new degenerate solutions of (|3~K1) appear, for both j = and j = 1, corresponding to 
the birth of a pair of fixed points, one stable and the other unstable. In Figure 2 we plot 
the left hand side of (|35|) with e = 2.2 and N = 1597, as above, but now with k = K*. It is 
apparent that the scars around the two bifurcations, at positions q = 0.05 and q = 0.44, 
are wider than those around the two non-bifurcating fixed points, at positions q = 0.69 
and q = 0.81, and that they are also wider than those shown in Figure 1. In Figure 3 we 
plot the left hand side of fpTHt) close to the bifurcation point at q = 0.44, together with the 
approximation (|39|) . which clearly captures the details of the associated fringe structure. 

It is straightforward now to deduce the semiclassical scaling with N of the moments 

C 2m (N)=N 2m J2 ((\*n(Q)\ 2 )-T7) ■ (41) 
Q=l ^ ' 

The arguments of Section 3 suggest that C<im is of the order of jV~ m+w , where w = 1/2 
away from bifurcations and w = 2/3 at the bifurcation (w is one plus the width exponents 
deduced from (f39|) . because those were for y rather than Ny, as we need here). In Figure 
4 we plot logC 2 against log N when k = 3. The fact that the points lie on a straight line 
confirms that there is a power-law scaling; furthermore, the gradient is close to —1/2, as 
expected. In Figure 5, we make the same plot for k = n*. In this case the gradient is close 
to the expected value of —1/3 (a possible explanation for the deviation is that for the 
range of values of iV shown, the bifurcation exponent is contaminated by the contributions 
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from the non-bifurcating periodic orbits). Finally, in Figure 6 we plot 



log C 2m {N) 

aim) = lim — — — , (42) 

yv ; n-,oo logN v ' 

calculated numerically from the gradients of best-fitting straight lines to plots like those 
in Figures 4 and 5. For both k = 3 and k = k* the results are in accord with the scaling 
law suggested above. 

We emphasize that these numerical computations illustrate the influence of one in- 
dividual bifurcation only. They do not test the competition which would result from 
averaging over a parameter range that contains many different generic bifurcations, and 
which the analysis of Section 3 suggests has a universal outcome for the moment expo- 
nents. 
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Figure 1: ^ n |^„(Q) | 2 5 £ (0 - n ) - 1, with e = 2.2, k = 3 and JV = 1597 (circles connected 
by dotted lines). Also shown is a convolution of the data with a normalized Gaussian of 
width 0.02 (bold line). The positions of the fixed points are q = 0.65 and q = 0.85. 
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Figure 2: £ n \ip n {Q)\ 2 4 (9 - 6 n ) - 1, with e = 2.2, k = k* and iV = 1597 (circles 
connected by dotted lines). Also shown is a convolution of the data with a normalized 
Gaussian of width 0.02 (bold line). There are unstable fixed points at q = 0.69 and 
q = 0.81, and bifurcations at q = 0.05 and q = 0.44. 
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Figure 3: J2 n \^n(Q)\ 2 S E (9 - 9 n ) - 1, with e = 2.2, k = k* and N = 1597, as in Figure 2, 
in the neighbourhood of the bifurcation at q = 0.44 (circles connected by dotted lines). 
Also shown is the local approximation obtained by substituting (|39|) into (|35|) (bold line). 
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Figure 4: logC2, calculated using a local Q-average of size 0.02iV 1/ ' 2 , plotted against logiV 
when k = 3 (circles). Also shown is a best-fitting straight line, which has gradient -0.50. 
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Figure 5: logC2, calculated using a local Q-average of size 0.02iV 1/ ' 2 , plotted against logiV 
when k = k* (circles). Also shown is a best-fitting straight line, which has gradient -0.29. 
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Figure 6: g(m) plotted against m, for k = 3 (full circles), and k — k* (open circles). Also 
shown are the best fitting straight lines: g = —0.96m + 0.57 and g = —1.01m + 0.72, 
respectively. 



24 



